Self-organized transient facilitated atomic transport in Pt/Al(lll) 



OO 

o 
o 

(N 



C/3 



P. Siile 

Research Institute for Technical Physics and Material Science, 
Konkoly Thege u. 29-33, Budapest, Hungary, 
sule@mfa.kfki.hu, www. mfa. kfki. hu/^sule, 

(Dated: February 5, 2008) 

During the course of atomic transport in a host material, impurity atoms need to surmount an 
energy barrier driven by thermodynamic bias or at ultra-low temperatures by quantum tunneling. 
In the present article we demonstrate using atomistic simulations that at ultra-low temperature 
transient inter-layer atomic transport is also possible without tunneling when the Pt/Al(lll) im- 
purity/host system self-organizes itself spontaneously into an intermixed configuration. No such 
extremely fast athermal concerted process has been reported before at ultra low temperatures. The 
outlined novel transient atomic exchange mechanism could be of general validity. We find that the 
source of ultra-low temperature heavy particle barrier crossing is intrinsic and no external bias is 
necessary for atomic intermixing and surface alloying in Pt/Al although the dynamic barrier height 
is few eV. The mechanism is driven by the local thermalization of the Al(lll) surface in a self- 
organized manner arranged spontaneously by the system without any external stimulus. The core 
of the short lived thermalized region reaches the local temperature of 1000 K (including few tens 
of Al atoms) while the average temperature of the simulation cell is 3 K. The transient facilitated 
intermixing process also takes place with repulsive impurity-host interaction potential leading to 
negative atomic mobility hence the atomic injection is largely independent of the strength of the 
impurity-surface interaction. We predict that similar exotic behaviour is possible in other materials 
as well. 

PACS numbers: 66.30. Jt, 68.35.Fx, 05.45.-a, 66.30.-h, 81.10.-h 
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I. INTRODUCTION 

The spontaneous formation of self-organized nanoscale 
structures has attracted much attention in recent years 
due to its potential application in the fabrication of 
nanodevicesii^. The understanding of atomistic pro- 
cesses which lead to the formation of nanostructures has 
been one of the main focuses of research in materials 
scienceii^. 

The self-assembly induced atomic movements towards 
an ordered structure can be understood as thermally acti- 
vated processes driven by the thermodynamic bias of the 
system^'**. Also, concerted atomic transport processes 
during self-organization such as adatom nucleation via 
detechment and attachment processes at step edges and 
thin film growth and processing, however, often lead to 
abrupt surface alloying and intermixin^i^iii^. These pro- 
cesses proceed via atomic site exchanges within the top- 
most atomic laye r^i^°i^^i^^ . 

Ultrafast diffusional dynamics can be studied by clas- 
sical molecular dynamics (MD) simulations at the atom- 
istic level— ii^ii^iiiii^ii^. Reccntty it has been shown by 
MD studies in accordance with experimental results that 
under externally forced conditions, transient enhanced 
intermixing of heavier impurities could occur in bulk 
materials In the absence of considerable external 

load of perturbation, such as during atomic deposition 
ultrafast intermixing and surface alloying can also be 
induced^'^'i^ii^'^i. Moreover, the most recently it has 
also been found that the bulk mobility of atomic metal- 
lic clusters could also be enhanced leading to ballistic 



burrowing in Al and in Ti2S. 

These are interesting results beacuse it is widely ac- 
cepted that enhanced diffusion occurs the mostly on 
solid surfaces e.g. when the barrier of atomic transport 
AE < ksT, where and T are the Boltzmann con- 
stant and the temperature, respectively, superdiffusion 
occurs, that is the nearly dissipationless atomic trans- 
port with transient atomic jumps (random walk. Levy 
flight In the topmost layer fast atomic ex- 
change processes with long jumps have also been re- 
ported which thought to be driven, however, by ther- 
modynamic forces—iiiii^. In the bulk, non-Arrhenius 
(athermal) atomic (not necessarily transient) transport 
has been studied mostly under nonequilibrium condi- 
tions, such as during ion-implantatioi>i^, in driven-alloys 
or mechanical alloyin g^^'^^ , by mechanical force biased 
chemical reactions^! or using shock- induced alloying'*^. 

In the bulk athermal rates can only be accom- 
plished via under barrier atomic jumps called quantum 
tunneling^i^^. This can be done mostly for light par- 
ticles at ultra-low temperatures. Quantum diffusion of 
H has been studied in detail on solid surfaces^". In 
the bulk, only few light elements show ultrafast inter- 
stitial diffusioni i^^'^^ . However, the quantum tunneling 
diffusion of heavy adatoms on various substrate surfaces 
have also been observed recentl y'^^i'^^ . Quantum diffusion 
(QD) have not been observed yet for heavy elements in 
the bulk although the de Broglie wavelength could be in 
the range of tunneling distance which allows QD on the 
surface^. The most recently reactive diffusion dynamics 
has been interpreted as a superdiffusive process during 
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the front propagation of interfaces^ which could be the 
first (though theoretical) finding that transient diffusive 
atomic transport takes place in the bulk. 

We would like to present classical MD results which 
suggest that transient rates could also be occurred with- 
out tunneling at ultra-low temperature via a pecu- 
liar mechanism during surface alloying. The employed 
semiempirical approach is validated by ab initio den- 
sity functional calculations. The explored new transient 
atomic exchange process is driven in a self-organized 
manner: the impurity/host system spontaneously reor- 
ganizes itself in such a way that abrupt surface alloying 
takes place at an ultra-low temperature. We also find 
that impurity induced local thermalization occurs at ~ 
K external temperature although no forced condition has 
been applied. Such a mechanism has not been observed 
yet although is likely to be of general validity. The local 
thermalization of the substrate facilitates atomic injec- 
tion while the overall temperature of the substrate re- 
mains very low (few K). A surprising consequence is that 
repulsive intermixing (and negative transient atomic mo- 
bility) could also occur theroretically, the process is not 
sensitive to the strength of the Al-Pt interaction. 



II. THE SIMULATION APPROACH 

Classical tight-binding molecular dynamics 
simulations^ were used to simulate soft landing 
and vapor deposition of Pt atoms on Al(lll) substrate 
at ~ K using the PARCAS codei? which has been used 
for the study of various atomic transport phenomena in 
the last few years-'^Q. We also employ first principles 
calculations to validate our heteronuclear potential 
(details will be given later on). 

Although we carry out simulations at ~ K, we find 
a substantial local heating up in a local surface region 
of Al, hence the correct dissipation of the emerged heat 
should be handled via using temperature control. A 
variable timestep and the Berendsen temperature con- 
trol is used at the cell border— i^i^. The simulation 
uses the Gear's predictor-corrector algorithm to calcu- 
late atomic trajectories^^. The maximum time step of 
0.05 fs is used during the operation of the multiple time 
st ep algorithm^. The system couples to a heat bath 
via the damping constant to maintain constant tempera- 
ture conditions and the thermal equilibrium of the entire 
system^. The time constant for temperature control is 
chosen to be r = 70 fs, where r is a characteristic relax- 
ation time to be adjusted^i^. The Berendsen temper- 
ature control has successfully been used for nonequilib- 
rium systems, such as occur during ion-bombardment of 
various material a-'^^'^^i-'^^'^'^i'^'-''^"'^ . Further details are given 
in ref42ii^ and details specific to the current system in 
recent communications^ii^i^. 

For simulating deposition it is appropriate to use tem- 
perature control at the cell borders. This is because it is 
physically correct that potential energy becomes kinetic 



energy on impact, i.e. heats the lattice. This heating 
should be allowed to dissipate naturally, which means 
temperature control should not be used at the impact 
point. Periodic boundary conditions are imposed later- 
arily. The observed anomalous transport processes are 
also observed without periodic boundary conditions and 
Berendsen temperature control. Further details are given 
in"*^ and details specific to the current system in recent 
communicationsi^i^. 

The top of the simulation cell is left free (the free sur- 
face) for the deposition of Pt atoms. The bottom layers 
are held fixed in order to avoid the rotation of the cell. 
Since the z direction is open, rotation could start around 
the z axis. The bottom layer fixation is also required to 
prevent the translation of the cell. 

The size of the simulation cell is 80 x 80 x 42 A'^ in- 
cluding 16128 atoms (with a fee lattice). The simulation 
uses the Gear's predictor-corrector algorithm to calcu- 
late atomic trajectories^^. The maximum time step of 
0.05 fs is used during the operation of the multiple time 
step algorythm"*^ . 15 active MLs are supported on 3 fixed 
bottom monolayers (MLs). We find no dependence of the 
anomalous atomic transport properties of the deposited 
atoms on the finite size of the simulation cell. Finite 
size effects do not play a role in the appearance of the 
anomalous transport of Pt in Al (the variation of the cell 
size does not influence the intermixing process down to 
cell sizes including few hundreds of atoms). Deposited 
atoms were initialized normal to the (111) surface with 
randomly selected lateral positions 4 — 5 A above the 
surface with nearly zero velocity. The initial kinetic en- 
ergy of the deposited particles in the case of ultrasoft 
landing is nearly zero eV. In order to make a statistics 
of impact events we generated 100 events with randomly 
varied impact positions. The conservation of the total 
energy is maintained during the simulations. 

A. The interaction potential 

We use the many-body tight-binding second-moment 
approximation (TB-SMA) interaction potential to de- 
scribe interatomic interactions'^^. Using the Cleri-Rosato 
(CR) parameterization of the TB-SMA potential we con- 
sider the interaction between two atoms and the interac- 
tion with their local environment. 

The TB-SMA potential is formally analogous to the 
embedded atomic method (EAM;^) formalism, e.g. the 
potential energy of an atom is given as a sum of repulsive 
pair potentials for the neighboring atoms (usually for the 
first or second neighbors and a cutoff is imposed out of 
this region) and an embedding energy that is a function 
of the local electron density given as follows^, 

Etot = lY.^{n,) + J2F[p,], (1) 

ij i 

where is the distance between atoms i and j. There 
are many functional forms are available for the density 
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FIG. 1: The crosspotential energy (eV) for the Al-Pt dimer as 
a function of the interatomic distance (A) obtained by the ab 
initio PBE/DFT method. For comparison the interpolated 
semiempirical potential (TB-SMA) is also shown calculated 
for the Al-Pt dimer. 



Pi and for the embedding function F[pi]^. In the code 
PARCAS^ the forces have been calculated using a built- 
in functional derivative of Eq. (1). We utilize EAM func- 
tional forms in the code for F[pi\ and for the density p 
similar to that given in 

j,gfg^45,53^ The EAM routine in the 
code employs a cubic spline interpolation for the evalua- 
tion of the EAM potentials and their derivatives (forces) 
starting from various kind of input potentials given in 
discrete points as a function of (the number of points 
per functions is 5000 in this study). 

Within the TB-SMA we use no explicit dependence on 
Pi. The attractive part of the potential reads, 



exp 



j-rij<ra 



-2q ( ^-1 



1/2 



, (2) 



where Tc is the cutoff radius of the interaction and tq is 
the first neighbor distance (atomic size parameter). 

The repulsive term is a Born-Mayer type phenomeno- 
logical core-repulsion term: 



exp 



-P 



ro 



(3) 



The parameters {^,q, A,p,ro) are fitted to experimental 
values of the cohesive energy, the lattice parameter, the 
bulk modulus and the elastic constants cn, Ci2 and 044^^ 
and which are given in Table 1. The summation over j is 
extended up to fifth neighbors for fee structures^. The 
cutoff radius Tc is taken as the third neighbor distance 
for all the interactions. We tested the Al-Al and the 
Al-Pt potential at cutoff radius with larger neighbor dis- 
tances and found no considerable change in the results. 
This type of a potential gives a very good description 
of lattice vacancies, including migration properties and 
a reasonable description of solid surfaces and melting22.. 
The CR potential correctly provides the adatom binding 
and dimerization energies^i. 



TABLE I: The Cleri-Rosato parameters^ used in the tight 
binding potential (TB-SMA) give n in Eqs . (l)-(2)^ The 
parameters of the crosspotential have been obtained as follows 
using an interpolation scheme^: For the preexponentials ^ 
and A we used the harmonic mean AAipt = {Aai x Apt)^^"^ 
(^ has been fitted to the heat of mixing of the AlPt alloy 
phase, see details in ref.— ), for q and p we use the geometrical 
averages: qAiPt = {qAi +qpt)/2. The first neighbor distance 
of the Al-Pt potential is given also as a geometrical mean of 
ro = {rP + r^')/2. 
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ro 


Al 


1.316 4.516 


0.122 


8.612 


2.87 


Pt 


2.695 4.004 


0.298 


10.612 


2.78 


Al-Pt 


2.7 3.258 


0.191 


9.612 


2.83 



Recently it has also been shown, that the CR potential 
remarkably well describes diffusion in liquid Al^>^ and 
energetic deposition of Al clusters on Al^. For the Al-Pt 
crosspotential of substrate atoms and Pt we employ an 
interpolation scheme which has widely been used in the 
literatureiSii^ii^iSii^i^I. The Al-Pt potential provides a 
reasonable melting point and heat of alloying for the AlPt 
alloy 

In order to check the accuracy of the employed inter- 
polated crosspotential, the crosspotential energy has also 
been calculated for the Al-Pt dimer using ab initio lo- 
cal spin density functional calculations-^ together with a 
quadratic convergence self-consistent field method. The 
G03 code is well suited for molecular calculations, hence 
it can be used for checking pair-potentials. The Kohn- 
Sham equations (based on density functional theory, 
DFT)^^ arc solved in an atom centered Gaussian basis 
set and the core electrons are described by effective core 
potentials (using the LANL2DZ basis set)^ and we used 
the Perwed-Burke-Ernzerhof (PBE) gradient corrected 
exchange-correlation potential^. First principles calcu- 
lations based on density functional theory (DFT) have 
been applied in various fields in the last few years^. 

The obtained profile is plotted in Fig. [1] together with 
our interpolated semiempirical many-body TB-SMA po- 
tential for the Al-Pt dimer. We find that our interpolated 
TB-SMA potential when calculated for the Al-Pt dimer 
matches reasonably well the ab initio one hence we are 
convinced that the TB-SMA model accurately describes 
the heteronuclear interaction in the Al-Pt dimer. We as- 
sume that this dimer potential is transferable for those 
cases when the Pt atom is embedded in Al. This can 
be done because, as we outlined above, the interpolated 
Al-Pt potential properly reproduces the available exper- 
imental results for the Al-Pt alloy. 

We also calculate the binding energy [/^ of the impurity 
particle i which can be expressed in terms of its potential 
energy U{r) — (summed over interactions with its neigh- 
bors cut off at Tc) and its first derivative (Newtonian 
forces). Hence at each time step UJ^ can be calculated 



FIG. 2: The snapshots of the uhrafast atomic injection at 
0, 0.45, 0.55, and at 2, ps (from left to right) 0.8, 1.2 and at 
1.5, ps (from left to right) Pt atom and Al atoms are shown 
with light and dark (blue) colors, respectively. Only few hun- 
dreds of atoms of the simulation cell in the "active" region 
are shown. The deposited particle is in rest at f = K (no 
initial kinetic energy is given). Therefore the atomic injection 
occurs spontaneously. 



from the Newtonian forces. 



E 



dUjr) 
dr 



(4) 



j,rij<rc,i^j ' 

Since the system is energy conservative, the space inte- 
gral over the Newtonian interatomic forces ^^^^ 

give the total energy of the simulation cell, r can be 
replaced by the internuclear separation r^j in the pair 
interaction term U{rij). 



III. RESULTS 

Transient inter-layer atomic mixing (TILAM): The 
simulation of vapor deposition in the Pt/Al(lll) system 
leads to an unexpected result. The deposited atoms, in- 
dependently of the energy of deposition, intermix sponta- 
neously with ultrafast atomic exchange entering the top 
Al(lll) layer even at '--^ K within a ps leading to the ex- 
tremely large jumping rate of F « 10^^ Hz. Such a robust 
rate at '--^ K has never been reported before in metal 
system. Moreover we find that the deposited Pt atom 
undergoes an abrupt inter-layer migration spontaneously 
regardless to the impact energy and to the strength of the 
Pt-Al interaction potential. Hence we find it important 
to understand the details of this porcess. Although the 
obtained rate is surprisingly high, the employed simula- 
tion approach is highly standard and hopefully there is no 
reason to question the validity of the results. In particu- 



FIG. 3; The vertical positions of moving atoms as a crossec- 
tional view of typical trajectories of Al atoms in the upper 
layers (in a crossectional slab cut the middle of the simula- 
tion eel 1) induced by the deposition of a Pt atom at ~ K 
(the scale is in A in the axes, the depth position of the surface 
is at z = 0). The trajectory of the Pt is also shown with a 
black curve. The positions of the atoms are collected up to 
2 ps during a deposition event. The different colors of the 
points correspond to the local temperature range (K) of the 
Al atoms in the thermalized region shown above. Inset on the 
top: The fluctuating local average temperature {Tiocai{t), K) 
of the thermalized region (~ 10 x 10 x 5 A^) which includes 
30 — 40 hot atoms as a function of the time (ps). 



lar, we do not think that the result is the artifact of the 
employed semicmpirical interaction potential. Our DFT 
dimer calculations support the reliability of our interpo- 
lated semicmpirical potential. We find such a peculiar 
behavior only for few diffusion couples (Pt/Cu, Pt/Al, 
Au/Al) among those couples for which interpolated CR 
potential has been available. Many experimental results 
support indirectly our finding. Strong exothermic solid 
state reactions have long been known between various 
metals and Al^'^^.so^ 

The injection of a Pt atom leads to the ejection of an Al 
atom to the surface. The deposition of 1 ML of Pt leads 
to the formation of an adlayer rich in Al in agreement 
with the experimental findings^i^^. The available room 
temperature experimental results also report us strong 
intermixing for Pt/Al^i^. The computer animation of 
the atomic injection can be seen in a web page^. We 
find direct injection only in the case of certain transition 
metal elements around Pt in the periodic table, such as 
Ir, Au and also for the Pt/Cu couple. 

The spontaneous local thermalization of the substrate: 
In order to get more insight into the details of the atom- 
istic mechanism of TILAM we follow the atomic trajec- 
tories of surface Al atoms. The TILAM induced surface 
disordering of Al(lll) is shown in Fig. [3] where the tra- 
jectories of the transient vertical jump of few Al atoms 
to the surface (adlayer) can be seen. We plot the atomic 
positions of a crossectional slab cut in the middle of the 
simulation cell (with a slab thickness of 15 A) for the 
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top layer atoms during vapor deposition of Pt (Fig [3]) 
at '--^ K. The mobility of the substrate atoms is large 
with large amplitudes around their equilibrium positions. 
This is surprising because the overall temperature of the 
entire cell does not exceed few K during the simulations. 
However, we find that within a small volume below the 
surface including few tens of atoms the local temperature 
can be surprisingly high. 



A. Local temperature within a thermalized 
subsurface nanoscale zone 

An approach will be outlined briefly which has been 
used for obtaining continuously distributed local proper- 
ties from the positions and velocities of constituent atoms 
obtained by MD simulations. The applied methodology 
is similar to that obtained for analyzing the results of MD 
simulations for systems with finite size and which can not 
be described by continuum models^i^. Unfortunatelly, 
when the system size is shrinked to the nanoscale fluc- 
tuations, such as the spatial oscillation of the tempera- 
ture will be enhanced^. Therefore specific definitions are 
needed for giving time dependent local quantities, such 
as the local temperature of a nanoscale system. 

Thermondynamic quantity, such as the temperature T 
can only be assigned to an atomic ensemble in which the 
number of particles N is sufficiently high to exclude the 
effect of local fluctuations (statistical ensemble average). 
Unfortunatelly, a nanosystem, which often includes less 
than ~ 1000 atoms (nanoclusters) can not be described 
by T in a conservative point of view. In these cases, how- 
ever, one should introduce the quantity local temperature 
Tiocai which can be used to explain the thermal properties 
of nanostructures. We define {Tiocaij as the time aver- 
aged temperature of the thermalized sub(nano)region ob- 
tained during simulations sampling the sufficiently large 
portion of the phase spaced. As a natural consequence, 
{Tiocai) T, when iV ^ oo, or is sufficiently large to 
have a statistical meaning of the ensemble. 

The total simulation cell of Pt/Al during the inter- 
mixing of Pt is a highly anisotropic and inhomogeneous 
system. In this case the thermalized region of the sub- 
strate can be taken as a nanosystem including few tens of 
hyperthermal atoms at the surface, however, which is not 
isolated from its low-temperature environment. There is 
a continous thermal exchange between the "hot spot" 
and the ultra-low temperature environment and with the 
heat bath. 

Much effort has been put forward in establishing 
a relationship between thermodynamic quantities and 
MD datapi^i^i^i^i^. The ergoditic theorem ensures 
in thermodynamics the relation between the observ- 
able ensemble averages and the simulated time averaged 
quantities^i^. Under ergodic condition the time average 
or ensemble average of the velocity distribution of the 
constituents will closely follow the Maxwell-Boltzmann 
distribution. MD studies provides microscopic (atomic) 
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FIG. 4: Ftg 4a: The local temperature Tiocai{t) (K) of the 
thermalized region below the surface (with the volume of 
10 X 10 X 5 A'^) as a function of time (ps) obtained for a 
typical event. The average temperature of the simulation cell 
is also shown with a dashed line. Inset: The number of AI 
atoms in this subsurface nanoscale region as a function of 
time (ps). Fig 4b: The oscillating kinetic energy (eV) of the 
impinging Pt atom and of the transient Al atoms (summed up 
for the thermalized ensemble) as a function of the time (ps) 
shown with a continous and dashed lines, respectively. The 
Pt atom has been initialized 4.6 A above the Al(lll) surface 
with zero velocity. The initial temperature of the Al cell is 
nearly zero. Fig 4c: The average cohesive energy/atom (eV) 
in the thermalized region of Al as a function of time (ps). 
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information and an effective temperature can be derived 
from the individual atomic velocitiesi^i^. Another use- 
ful quantity, the kinetic temperature of a nanosystem is 
defined as an averaged kinetic energy per an individual 
spatial degree of freedom^^'^'^. The time averaged tem- 
perature of the nanosystem with N number of atoms can 
be given as 

1 /"^ 

{Tlocal)^ = lim - / Tiocal{t)dt 
t Jt=0 

= lim - / — '^^P^dt 

t Jt=0 3fcB 

where M is the number of visited configurations in the 
simulations. The instanteneous value of Tiocai (t) will fluc- 
tuate around the mean value {Tiocai)^ unless an infinite 
number of particles haven been considered. Indeed, un- 
der the assumption of crgodicity, this mean value exists 
and is independent of the initial data with identical en- 
ergy. In nanosystems ergodicity looses its validity and we 
can not give any macroscopic observable which is related 
to (Tiocai)- Nevertheless, the time evolution of Tiocai(t) 
could be a useful quantity to monitor system changes dur- 
ing transient phase evolutions such as e.g. local melting 
transition'^^. 

The calculation of Tiocai (t) within an arbitrarily small 
volume including N particles can also be calculated for- 
mally, however, from the kinetic energy of the individual 
particles obtained from simulations in each time step. 
The equipartition theorem (ET) allows in principle to re- 
late the temperature of a system with its average energy. 
If the mobility of atoms is sufficiently high within the 
thermalized nanoregion, e.g. the number of hyperther- 
mal atoms is large, we are close to the limit of classical 
ideal gas. In this case each of the mobile particles has an 
average kinetic energy of (3/2)fcsr in thermal equilib- 
rium, where ks is the Boltzmann constant and T is the 
temperature. However, our system should not be strictly 
an ideal gas. ET requires only that k{Tiocai) 3> hv, where 
f is the phonon frequency of an oscillator in the thermal- 
ized zone, hence in highly excited states quantum effects 
should become negligibleSi. 

Although our system is not in a thermal equilibrium, 
we can expect that the ET works also nearly correctly 
for those cases which are not very far from equilib- 
rium. It is also known that the spatial-temporal vari- 
ation of Tiocaiit,r) can be calculated for even highly 
inhomogeneous systems such as e.g. a plasma, shock 
loaded systems^ or for ion-bombardment induced ther- 
mal spikes^^. The thermalized region if assumed as an 
overheated liquid state of matter, and is close to an ideal 
gas and than ET can be applied. In an ideal gas atoms 
can move few A/ps (that is a nearly ballistic atomic mo- 
bility). In our system we find a similar rate of mobility 
within the thermalized region during simulations when 



a Pt impurity atom has been injected. The thermal- 
ized state of the local subregion in Al persists up to 
few ps which is a very short lifetime in the thermody- 
namic sense. The quench rate of the liquid-like phase 
is extremely fast. Nevertheless, the local temperature 
Tijocai{t) can be formally assigned to each particles dur- 
ing the molten phase at time t using the ET, 

1 3 

-mivf{t) = -ksTijocaiit), (6) 

where Vi{t) and nii are the atomic velocity and mass of 
ith particle. Note that we do not follow the spatial vari- 
ation of Tijocai{t) as it has been calculated e.g. in refi^. 
Ti^iocai{t) is always assigned to the velocity Vi{t) of par- 
ticle i at time t. We are interested in the time evolution 
of Ti^iocaiit) in a particle trajectory. Ti^iocai{t) has no 
physical meaning in a strict thermodynamic sense since 
ergodic theory relates only the time averaged (Tiocai) to 
the ensemble average (observable temperature). How- 
ever, summing up for each particles within the thermal- 
ized volume and calculating the average Tiocai{t) of an 
ensemble of particles in a local volume could provide 
a time dependent ensemble averaged local temperature 
with physical meaning. The averaged local temperature 
Tiocai (t) within a region of the substrate is given than at 
time t for N number of hyperthermal atoms 

Tlocal{t) = ^ X! ^^^localit) = ^ X! ' 

We use Tiocai (t) for describing the time evolution of local 
heating up (thermalization) processes. Probing few tens 
of configurations with different starting positions for Pt 
we get very similar events hence we do not average for 
events in Eq. [T] Hence the plotted Tiocai{t) values do 
not correspond to an ensemble average, instead we show 
the time evolution of the simulated Tiocai{t) for a typi- 
cal Pt/Al deposition event in Fig. 3^. Therefore the 
plotted Tiocai{t) curve is not a unique one, the shape of 
temperature fluctuation changes event by event. How- 
ever, in most of the events we find similar features: mul- 
tiply peaked structure (transients) and high average local 
temperature within the subsurface zone {Tiocai) ~ 500 K 
during the deposition. Such kind of an analysis can be 
used e.g. for studying the time evolution of local melting 
in coUisional cascades and thermal spikesi^i^^. 

The magnitude of Tiocai (t) clearly depends on the vol- 
ume of the subregion considered in the summation in 
Eq. [71 We choose a subregion in which the individual ki- 
netic energy of particles exceeds a threshold value (that is 
Ekin > ^kTraeit, whcrc Tmeit IS the bulk melting point of 
e.g. Al, E'^-'^ K. 0.07 eV/atom in this case). This thresh- 
old value is rationalized by our experience reported in 
recent publications^ in which we found that local melt- 
ing transition can be described by the occurrence of suf- 
ficiently large number of " liquid" particles which possess 
Ekin > ^kTraeit- Out of the subrcgiou we cut off con- 
tributions to Eq. [T] The crossection of the thermalized 
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volume at and below the surface can be seen in Fig. [3] 
which shows us that the volume of this region is less than 

10 X 10 X 5 A'^ including less than ~ 200 atoms. 

Using the local temperature analysis outlined above 
we find that in the core of the disordered region (20 — 40 

atoms in the subsurface zone of ~ 5 x 5 x 2.5 A ) the peak 
local temperature (LT) reaches Tiocai{t) ~ 1000 K for few 
times during the persistense of thermalization (see Fig. 

Hi). In the large zone of - 10x10x5 A^ Tiocai{t) ^ 100- 
150 K. This is really surprising that Tiocai can reach so 
high instanteneous values while the external temperature 
is ~ K. This is because no external forced condition 
has been applied: simply the presence of an impurity 
particle induces the local thermalization. This can be 
taken as a spontaneous process: no driven conditions has 
been applied. The impurity particle has been added to 
the host system in a rest since no initial velocity has 
been given to it. We are not aware of other reports in 
which the spontaneous local heating up of a K system 
has been published before. This finding could explain 
the strong exotermicity known during alloying between 
various transition metals and Al at room temperature 
which could even lead to extremely fast burn rates^. 




depth 

FIG. 5: The dynamical potential energy {{U[z)), binding 
energy) of the impinging Pt atom (eV) during a single atomic 
jump through the topmost layer of Al(lll) as a function of 
the distance from the surface (in A). The crossectional views 
of the thermalized region of the system are also shown with 
the impinging impurity particle for the initial system and for 
the transition (activated) structures (which correspond to the 
barriers) at t = 0, 0.85, 0.95, 1.05 and 1.15 ps. 



B. The driving force of TILAM: kinetic energy 
transfers 

We also find that the system exhibits a complicated 
time evolution. We plot various local properties against 
time, such as Tiocai{t), the number of Al atoms in the 
thermalized subregion (Fig. |4^), the kinetic energy of 
the Pt atom and the Al atoms (the sum of the kinetic 
energy of the thermalized Al atoms. Fig. |3l3) and the 
average cohesive energy in the thermalized nanoregion 
(Fig. Ht). 

We notice in Figs. |4] the multiply peaked structure 
of the curves. The kinetic energy of Pt oscillates as a 
function of time and which correlates with the fluctua- 
tion of the LT of the thermalized region of Al shown in 
Fig. |4^. The correlation also holds with the average 
kinetic energy of the hot Al atoms shown in Fig. In 
particular, we see 4 nearly equidistant peaks separated 
by ~ 0.1 ps in Fig. [3]d. This peculiar feature of the 
system does not support a heat spike mechanism (that is 
the simple impact induced thermalization and local melt- 
ing). The low impact energy of ~ 2 eV alone (obtained 
during local acceleration) is insufhcient for the appear- 
ance of a collisional cascade not even for a thermal spike 
(the formation energy of a vacancy-displaced atom pair, 
Frenkel pair is ~ 25 eV). There must be an additional 
mechanism which thermalizes the top region of the Al 
substrate. 

We would like to emphasize the spontaneous nature 
of the local heating up because we let the system to 
time evolve during the simulations without any external 
perturbation. No kinetic energy is given to the impu- 



rity atom, hence no externally forced condition has been 
applied. However, the system reorganizes itself sponta- 
nenously. This is typically the characteristics of a self- 
organizing system. 



C. The dynamic adsorbate-surface potential energy 
profile: 

In Fig. O we show the dynamic potential energy pro- 
file of the approaching Pt atom together with the crossec- 
tional views of the structures corresponding to the energy 
barriers of the oscillatory potential. The potential energy 
of Pt as a function of the distance from the surface {U (z)) 
has been calculated using Eq. 31 The barrier height of 
these peaks is around few eV. 

The static Pt-Al(lll) potential energy profile (and sur- 
face), which could be calculated e.g. by ab initio den- 
sity functional approaches, does not account for the dy- 
namic nature of TILAM. We see from the animations— 
the strong mobility of the substrate atoms during the ul- 
trafast atomic exchange. Hence the dynamic Pt-Al(lll) 
potential energy profile, obtained by MD, describes more 
correctly the energetics of the system. 

In Fig [5] we show the binding energy (potential en- 
ergy) of the Pt atom as a function of the distance from 
the surface. The binding energy has been calculated as 
a sum of pair interactions between the Pt and Al atoms 
up to the second neighbors. Further neighbor interac- 
tions (3rd and higher orders) have vanisingly small con- 
tribution. For typical events multiply peaked oscillatory 
impurity-substrate interaction potential have been found. 



8 




temperature (K) 

FIG. 6: The thermally activated jumping frequency (m^/s) 
as a function of Tiocai obtained by using the expression of 
r = roea;p(— ■^) at different activation energies given in eV. 



It must be emphasized that the profile of the potential 
is not unique, e.g. the depth positions and the height of 
the barriers slightly scatter in various events, neverthe- 
less, the oscillatory behavior typically occurs for all of 
the events during which injection take place. 

The magnitude of the scatter does not exceed few 
tenths of eV for the barrier height and few tenths of A for 
the depth positions of them. Due to the nonuniqueness of 
the oscillatory behavior it is no use to plot a profile over 
a statistical average of various events. It must also be 
noted that a potential energy profile for a normal depo- 
sition event without injection (when the impurity atom 
becomes an adatom) is a unique single well potential. 

Oscillatory interaction profiles have already been re- 
ported on the surface of metals between adatoms^^ , how- 
ever, no reports have been found for IM. The deposited 
particle could go through the multiply barriered potential 
due to the kinetic energy obtained during the local accel- 
eration towards the surface and due to a still unknown 
mechanism which will be characterized in the next sec- 
tion. The magnitude of the peak kinetic energy (arrival 
energy) is in the range of few eVs. 

The particle can gain the adsorption energy (conden- 
sation energy) as a kinetic energy during the accomo- 
dation process on the surface. Sufficiently high kinetic 
energy gain allows the transient lateral movement of the 
particle on the surface until sticking (trapping) occurs-. 
In Pt/Al, however, we find that instead of the lateral 
transient mobility the deposition of the impurity particle 
leads to transient IL atomic mobility. 

IV. DISCUSSION 

A. Quantum tunneling diffusion is ruled out: 

Since low-temperature and ultrafast atomic transport 
has been known as quantum tunneling it could also be 



that the impurity injection of Pt to Al(lll) is also QD. 
In principle, classical molecular dynamics can account for 
QD implicitly via the adjusted the parameters, although 
no direct quantum effects has been incorporated into the 
interaction potentials. E.g. the QD of H in various met- 
als can be simulated by classical MDf^. Therefore we 
check this process for QD. The de Broglie wavelength of 
a Pt atom A = 2nh/ \/3mkT « 1.7 A (during the sim- 
ulations we find T w 3 K in the simulation cell) which 
is somewhat less than the inter-layer hopping distance of 
~ 2 — 3 A. However, if we estimate the tunneling jump 
rate F, which is approximated by^^i^ 

where w = ^j2E^JmJP- , b is the barrier width (& « 1 
A). We get a vanisingly small Ft„„ « 10^^^ Hz, which is 
very far from our finding of Tmd ~ 10^^ — 10^'^ Hz (~ 1 
inter-layer jump/ps). We conclude that atomic transport 
with quantum tunneling can not explain the occurrence 
of transient inter-layer atomic mobility and there must be 
a peculiar mechanism which promotes impurity particle 
acceleration through the top layer of Al at ~ K. Using 
the Van't Hoff-Arrhenius expression (which can be used 
for classical atomic transport) of F = TQexp{—^) gives 
also zero thermal rate of F « Hz (the preexponential 
Fo w 10" Hz-*). 

The dynamic activation energy is taken from Fig. 4, 
from the dynamic potential energy profile of Pt along the 
reaction coordinate {Ea = 1 — 2 eV). However, if we take 
into account the local heating up {Tiocai ~ 1000 K of 
the region where the impinging Pt atom is intermixed) 
we get F K, IQ^ — 10® Hz which is also far from our MD 
simulations. Although in the hot core of the thermalized 
region we find few thousands of K temperature, how- 
ever this does not allow such an increase of the thermal 
jumping frequency. This has been demonstrated in Fig. 
[6l in which we plot the thermally activated T {Tiocai) as 
a function of the local temperature in the thermalized 
nanoregion. We show the obtained curves with 3 differ- 
ent activation energy values. It can be seen that even at 1 
eV activation energy, which is below our calculated value 
{Ea w 2 — 3 eV) , we get F much below the simulated 
Pmu ~ 10^^ — lO^'^ Hz in the reasonable temperature 
regime (T < 2000 K). 

The rate can only be in accordance with the simulated 
value within the Arrhenius picture if the prefactor Fq 
is increased by few orders of magnitude or at very high 
temperature of T « 10^ K. However, recent results rule 
out the occurrence of anomalous preexponential factorsi^. 
The serious deviation of the jump rate from the Van't 
Hoff-Arrhenius equation is highly unusual'* and has only 
been found for such cases when Ea ~ kT- which often 
lead to superdiffusion. The discrepancy between the ther- 
mal and transient rates could be accounted for assuming 
an auxiliary mechanism which facilitates the amplifica- 
tion of intermixing. In the rest of the paper we outline 
the details of such a possible mechanism. 
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FIG. 7: The snapshot of the transient facihtated intermixing 
(also called as TIL AM) process driven by the assistance of 
few surface Al atoms (which move out-of-plane coherently) 
at 0.95 ps. The displaced Al atoms are shown with a darker 
color. The impinging impurity particle (Pt) is shown with 
a lighter color on top of the displaced Al atoms. Only the 
central region of the top layer is shown. 




FIG. 8: The schematic view of transient facilitated inter- 
mixing in the coexisting solid-liquid like system. TR and UC 
denote the thermalized subsurface nanoregion {Tiocai ~ 1000 
K) and the ultracold substrate (few K). Opened arrows with 
dashed line show the trajectory of out-of-plane Al atoms 
which give downward momentum to the approaching Pt atom 
pushing down it below the top layer of Al(lll). 



B. The kinetic energy oscillation of Pt: 

The downward momentum transfer provided by the 
transient out-of-plane Al atoms on the surface surpasses 
the short-ranged impurity-host repulsion which promotes 
TILAM through an energeticaly unfavorable transtion 
state. The coincidence in time of peaks of the time evo- 
lution of the kinetic energy of the impinging Pt atom in 
Fig. IUd and the local temperature Tiocai{t) in Fig. [4^ in 
the thermalized region supports this explanation. More- 
over, we see in Fig that the sum of the kinetic energy 
of the transient Al atoms reaches its maximum with some 
time delay when compared with the kinetic energy of the 
Pt atom Fig. 0)3. The delay occurs 4 times indicating 
that the momentum transfer really occurs. 

First the approaching Pt particle gains kinetic energy 
during surface acceleration and gives the obtained en- 
ergy to the Al surface which in turn is heated up locally 
(a peak occurs with some time delay for the Al atoms). 
Than a second lower peak appears for the Al atoms to- 
gether with a higher peak for the Pt atom which supports 
the idea of momentum transfer between the transient Al 
atoms and the Pt atom. This process is repeated for few 
times periodically. 

In Fig [7] we show the appearance of few transient Al 
atoms on the surface ait k, 0.95 ps which catalyse the in- 
jection of Pt. In particular, 4 surface atoms move out-of- 
plane coherently and opens up a channel for TILAM. The 
average distance of these atoms from the surface is not 
larger than ~ 1 A. At a critical proximity of the Pt atom 
to the surface {dptAi < 3 A) strong, mostly out-of-plane 
atomic displacements of few Al atoms at the surface set 
in which leads to the injection of the particle. Without 
the out-of-plane surface instability no injection occurs. 



the particle becomes an adatom. Hence the concerted 
motion of few surface atoms in a self-organized manner 
opens up a channel spontaneously for injection. The tem- 
porally ejected and returning transient Al adatoms give a 
downwards momentum to the approaching impurity par- 
ticle: this kind of a kinetic energy transfer is so effective 
that the impurity particle can move through the top layer 
although the dynamic barrier height of this process is few 
eV. 



C. The microscopic mechanism of TILAM: 

The obtained results in this article might support the 
following mechanism: The impurity Pt atom approaches 
the surface of Al(lll) by surface acceleration. The de- 
posited particle is injected to the substrate by few out- 
of-plane vibrations of few surface Al atoms in 4 or more 
steps. The Pt atom gains the kinetic energy of the tran- 
sient Al atoms which promotes atomic intermixing. The 
schematic details of the process can be seen in Fig. |8l 
The schem shows us that few transient Al atoms at the 
surface move out-of-plane and when returning back give 
downward momentum to the approaching Pt atom in- 
jecting it below the surface. 

In Fig. [5] we see, that the energy difference between 
the initial and final configurations (between the energy 
minima) is ~ —0.5 eV which is also the driving force of 
TILAM (although it is insufficient alone for barrier cross- 
ing). Also the TILAM of Pt is driven not purely by the 
surface acceleration of Pt: first the impurity atom accel- 
erates towards the surface, than slows down and transfers 
its kinetic energy to few of the Al atoms (local heating 
up) and in the next period gains kinetic energy again 
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provided by the vibrating Al atoms and this is repeated 
for few times. Hence the ultra-low temperature IM of Pt 
is accomplished via a peculiar mechanism which includes 
mutual kinetic energy transfers between few substrate 
atoms and the impurity atom. The approaching Pt atom 
induces local heating up (thermalization) and in turn the 
transient Al atoms provide back kinetic energy to the 
Pt atom. This mutual transfer of kinetic energy takes 
place few times which amplifies inter-layer atomic trans- 
port substantially. The spontaneous interplay of kinetic 
energy transfers between the impurity and hyperthemal 
host atoms is a self-organized process: no external load 
of energy or perturbation is needed to initialize transient 
facilitated intermixing (TFI). TFI proceeds via the coher- 
ent movement of few atoms in the thermalized nanoscale 
subsurface zone induced by the impinging Pt atom. 



D. The possible role of atomic mass anisotropy 

The importance of mutual kinetic energy transfers dur- 
ing TFI (TILAM) is further evidenced by the simulation 
result that if we interchange atomic masses between Pt 
and Al (this can be done in MD without changing other 
parameters) and further increase the atomic mass of Al 
to mhost > 800 g/mol TFI can be suppressed. Surpis- 
ingly this very small atomic mass anisotropy of (5 < 0.03, 
{5 = mimp/rrihost, where rriimp and ruhost are the atomic 
masses of the impurity and host atoms, respectively) 
is required to stop intermixing of Pt while the natural 
5 « 7.5. Hence in Pt/Al we find that TFI is rather in- 
sensitive to the variation of the mass anisotropy. 

However, TFI can also be induced in other impu- 
rity/host couples with sufficiently large atomic mass 
anisotropy where rriimp S> rrihost- Indeed, if we set in 
artificially large (5 = 6 — 8 e.g. in Ni/Al or in Cu/Al 
which do not show TILAM with natural 5, we also ob- 
serve TFI. No TILAM occurs in other mass anisotropic 
systems, such as e.g. Au/Ni {5 ~ 3.4) or in Pt/Ti (5 ~ 4) 
at natural mass anisotropy. However, setting in (5 > 7 
in Pt/Ti, TILAM can also be induced. In Au/Ni we 
could not induce TILAM even with extremely large 5. 
Nevertheless, we conclude from this that 5 plays some 
role, however, not only 5 governs TFI as the prototypi- 
cal Pt/Al system shows it. It should also be noted that 
the role of 5 in various intermixing processes has already 
been studied in detail recently^-ii^i2ii2& and strong mass 
effect has been found during various ion-bombardment 
induced intermixing processes. 



E. Repulsive intermixing and negative mobility: 

One of the most intriguing consequence of the peculiar 
mechanism of TILAM is that since the process is purely 
kineticaly driven the strength of the Al-Pt interaction 
does not affect IM. We find that even if repulsive crosspo- 
tential is used, TILAM occurs. Tuning the strength of 



Al-Pt interaction the speed and the frequency of atomic 
injection is unaltered. Only very strong repulsive po- 
tential suppresses TILAM. The transient mobility of the 
impinging particle against the repulsive potential of the 
surface can also be understood as a transient negative 
mobility. This is similar to that reported by R. Eichhorn 
et al. for Brownian motion^* . Negative mobility of a par- 
ticle occurs when the response of a system to an external 
load (apphed bias) is opposite to the direction of this 
applied force. In our case if we place an impurity parti- 
cle (that is Pt atom) in rest above the repulsive surface 
of Al(lll) the particle begins to accelerate towards the 
surface, against the repulsive field providing a negative 
response and mobility. 

Hence in the facilitated transport system of Pt/Al(lll) 
the facilitated atomic transport is largely independent of 
the forces of interaction between the intermixing impu- 
rity atom and the substrate (host) surface. A conven- 
tional thermally activated inter-layer transport proceeds 
at the expense of the binding force developed in the tran- 
sition state. 



F. Possible indirect experimental evidences of 
TILAM 

(i) The occurrence of an Al adlayer upon monolayer 
Pt depositions^ which is in accordance with our finding 
that after the injection of Pt an Al atoms is released to 
the surface as an adatom. 

(ii) The strong exothermicity^ and large negative heat 
of mixing^ observed in various metal- Al reaction during 
thin film growth could be the macroscopic fingerprint of 
local thermalization observed in the nanoscale. 

(iii) The large intermixing length reported recently* also 
supports our findings: repeating for many times Pt de- 
position we also get a strongly intermixed film. 

(iv) The strong ion-bombardment induced intermixing 
in Al/Pt bilayer— is reproduced by our simulations^. 
These findings are reproduced at least qualitatively by 
our Al-Pt potential hence we are convinced that the het- 
eronuclear potential adequatively describes interatomic 
interaction between Al and Pt. 

Finally, the direct experimental confirmation of this 
new atomic transport mechanism could be carried out 
by ultra-low temperature scanning tunneling microscopy 
measurements for subnionolayer deposition events. 

V. CONCLUSIONS 

We predict an unprecendented rich dynamic phe- 
nomenon that has not been previously anticipated. We 
found a unique situation in which an impurity particle 
can move spontaneously through the barrier without re- 
flection and external energy income even at ~ K. 

The microscopic mechanism of the anomalously fast 
inter-layer transport of the deposited impurity atom at 
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'--^ K in few film/substrate couples (X/A1(111), where 
X=Pt, Au, Ir and in Pt/Cu( 111)) has been explained in 
detail. According to our knowledge, no ^ K transient 
heavy atomic surface alloying and intermixing has been 
reported before. We find that the Al surface behaves 
like an ultra- low temperature " atomic trap" . This kind 
of a mechanism could be general, although we find until 
now only few systems (impurity-substrate couple) which 
show transient inter-layer atomic mixing (TILAM) such 
as Pt/Cu besides Pt/Al (also Ir/Al and Au/Al). Also, 
no such superdiffusive behavior has been reported yet for 
intermixing. 

We find that although the mechanism of transient 
inter-layer atomic mobility is classical, the impinging im- 
purity particle moves through the potential barrier ki- 
netically, however, the transition rate does not depend 
on the external temperature, hence the process is ather- 
mal. The strong deviation from the Arrhenius law and 
the multiply barriered potential energy profile also im- 
plies an unconventional mechanism. 

The understanding of the mechanism could help to 
explain and classify various athermal atomic transport 
processes known in the literature such as superdiffu- 
sion on solid surfaces^ii^ii^i^, quantum diffusionii^ i^^i^^ , 
transient enhanced diffusion and intermixingi^i^, tran- 
sient cluster burrowing^^, ultrafast adatom island 
nucleation^'^, coherent ballistic displacement of atom&iS,, 
the cooperative enhancement of surface roughening^, 
the facilitated intermixing near step edges-^ or the su- 
perdiffusive mixing in plastically deformed solids^^. 

During few of these processes no direct external forced 
conditions haven been applied, or the athermal pro- 
cess occurs beyond the spatial range of the external 
stimulus^-iiii^iii^i^. One can classify them as self- 
organized facilitated transient atomic transport processes 
which lead to superdiffusion due to still unknown or not 
clearly established reasons. Until now no considerable 
effort has been made to understand the driving force 
of such athermal transient atomic transport processes. 
This piece of a work could help in resolving the mecha- 
nism of unconventional atomic transport which could be 
widespread in nature. 

The deposited particle gains kinetic energy during lo- 



cal acceleration (its initial kinetic energy is zero) towards 
the surface of the substrate and the arrival energy of few 
eV is eligible to overcome the barrier height of the inter- 
layer transport of interdiffusion. The impinging Pt atom 
induces the local heating up (melting) of the Al substrate 
(heat spike) and the occurrence of few transient Al atoms. 
The impinging impurity atom is driven through a mul- 
tiply barriered reaction pathway. The process is purely 
kinetically governed: transient uppermost layer Al atoms 
give downward momentum to the approaching Pt atoms 
which in turn is injected to the Al bulk. 

The mechanism of TILAM is similar to that of the 
facilitated passive transport of molecules through a lipid 
cell membrane. In both cases the process does not require 
extra energy in contrast to other passive transport pro- 
cesses in which thermodynamic bias or concentration gra- 
dient drives diffusion. During facilitated transport "car- 
riers" at the interface "catalyses" diffusion. In our case 
few hyperthermal Al atoms promote the process which 
move coherently out-of-plane and kicking down the Pt 
atom when returning back to their equilibrium position. 
During the course of TILAM, 4 transient Al atoms facil- 
itate the channeling of the Pt atom through the topmost 
layer of Al(lll). TILAM is possible due to the delicate 
interplay between the kinetic and vibrational degrees of 
freedom of the impurity/host couple which allows tran- 
sient inter-layer atomic transport at ~ K. 

These results together with a previous report^^ indi- 
cate that enhanced heavy particle (including heavy atoms 
and atomic clusters) bulk mobility could occur in few 
impurity/substrate couples due to an unprecentended 
mechanism. Surprisingly, transient negative atomic mo- 
bility is also possible with TILAM: the impinging particle 
can move through the top layer of Al(lll) against a re- 
pulsive surface potential. The exotic behavior of various 
impurity/host couples is shown to be largely dependent 
on the atomic mass anisotropy of the system while insen- 
sitive to the strength of the heteronuclear interaction. 
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